Radiative corrections to scalar masses in de Sitter space 
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We compute the radiative corrections to the mass of a test boson field in an inflating space- 
time. The calculations are carried out in case of a boson part of a supersymmetric chiral multiplet. 
We show that its mass is preserved up to logarithmic divergences both in ultraviolet and infrared 
' domains. Consequences of these results for inflationary models are discussed. 
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Quantum field theory in an inflating universe is thought to be the playground of the physical processes that took 
] place during the early Universe . However little has yet been explicitly computed for self interacting fields 2] . Any 
light bosonic fields with non minimal coupling that exist during an inflationary phase is bound to produce, together 

^ ^ with the inflaton field, superhorizon fluctuations that eventually might be visible in mechanisms such as isocurvature 

J> . modes generation from multi-field inflation curvaton models 0,13 1 bent trajectory inflationary models 

\^ ' 0,01 or modulated fluctuations 'Tl'|. If such fields are self interacting, for instance with a quartic potential, such 
OO , fluctuations might develop significant non-Gaussian features that in turn could be detected (see ref . 12,] for a review) . 
This would be only possible however if the radiative corrections do not render the particle too massive, e.g. with a mass 
' larger than the Hubble constant H, to develop any fluctuations at all. One question then raised by phenomenological 
investigations of this physics is whether the mass of a self interacting bosonic field can be protected against radiative 
corrections. 

. Of course one expects such a theory to be renormalizable but it implies that the fundamental bare theory has to be 
' fine-tuned. The radiative correction to the mass of a self interacting boson should indeed be naively of the order of 
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I. MOTIVATION 



(5m^ ^ AMpj (A being the coupling constant). This question has been investigated in 22] in the case of the inflaton 



field. In this case A is generically small since, for the inflaton field A Afpj has to be of the order of H^. In the case 
(-H we are interested in, however, we have no such constraints on A which can then be close to unity. For such values of 
the coupling constant the scalar field becomes too heavy to develop significant fluctuations unless the cancellation is 
precise over ten orders of magnitude - since H^/Mp^ ~ 10"^" during inflation. 
. On the other hand light scalar fields are expected to be associated with fermions as parts of super-multiplets in 
i—{ ' supersymmetric theories - such as in D or F-term inflation models [isL IT^. From the Minkowskian behavior, one 
5^ , expects the largest divergences coming from the fermionic and bosonic loops to cancel out. The aim of this paper is 
precisely to compute the radiative corrections to a test boson mass embedded in a chiral multiplet and when it lives in 
an inflating universe. The test model will be the Wess-Zumino Lagrangian in an expanding universe, not necessarily 
assumed to be de Sitter. The background evolution of the Universe is assumed to be driven by an other sector of the 
theory. 

We set up the formalism in section II. In section HI we compute the two-point function to get the one loop effective 
masses of the fields in any spacially flat Friedman-Lemaitre- Robertson- Walker (FLRW) spacetime. The ultraviolet and 
infrared behaviors are discussed in section IV. We show in particular that the infrared divergences can be apprehended 
in a classical approach. In the ultraviolet domain, the correction to the mass is found to be logarithmically divergent 
and proportional to the spacetime curvature and the coupling constant. It shows that masses of light scalar fields do 
not get a large contribution compared to the Hubble scale. Our results are put into perspective in the conclusion. 
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II. GENERAL SETUP 
A. A toy model 

We aim at studying radiative mass corrections in a generic Friedmann-Lemaitre- Robertson- Walker background with 
flat spatial sections where the background is assumed to be driven by another sector of the theory. We will focuss on 
the following Lagrangian 

^ = ^5^V*V^,^ + [i^V^ ^{x)\/^ -m]^- A^(<^i - *7'02)V^ - m \^\' - mX |^|' + v^*) - \^\'' (1) 

where (/? = l/\/2(</'i + i02) is a complex scalar field, t/" is a four-dimensional Majorana spinor and the 7" are the flat 
Dirac matrices. We also defined 

The kinetic term for fcrmions involves the Vierbein which are defined as 

where the coordinates C^{x) correspond to a locally inertial frame while are the non-inertial coordinates. One can 
check that the Vierbein behave like vectors under Lorcnz transformation (index a) and under coordinate transforma- 
tion (index /i). The metric g^u is related to the Minkowkian one rjap through 

9^.u{x) = V'' ^{x)Vf\{x)i^c.i3- (4) 

To obtain an action which is invariant under Lorentz group and coordinates change, the derivatives of a spinor have 
to be modified in the following way 

V„ = K.^(a;)V^ = y„^(a;)(5^-fr^) (5) 



where the spin connexion is defined by 



TAx)^]^a^PV^''{x)d^Vp,{x), (6) 
a°'<^ being the spinorial representation of the Lorentz group generators given in terms of the Dirac matrices 

a-P='-[r,l% (7) 

Since we will focus on spacially flat FLRW spacetimes, the metric should be written in conformal time such as 
<?/jiy — The Vierbein are simply expressed as V"' ^{x) = aS"^ ^ and defining the conformal fields as 

Lp — aLp (8) 

i) = a^/^^ (9) 

the Lagrangian takes the form 



(10) 

Recalling that the scalar curvature is given by i? = %a!' jc? ^ the bosons formally acquire an effective time-dependent 
negative mass square proportional to the curvature 

m2f, = -ii?<0 (11) 

while the fermions stay massless. As we shall see in the following, quantum corrections to the mass in the context of 
a cut-off regularisation lead to logarithmic divergences that are proportional to the curvature, as we should expect 
from ((TT|) . 

One can note that this Lagrangian is deduced from a Wess-Zumino model with the superpotential 

= 1^2^ + Ia^^ (12) 
^ o 

However, this is only a Wess-Zumino inspired toy model. Supersymmetric models in curved spacetime should be 
carefully dealt with since the parameter of a supersymmetric transformation should be space and time dependent 
|l5lll6j |. The background structure may forbid any coherent global supersymmetric approach. 
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B. Quantizations of the fields 



As we are interested in quantum corrections to masses, we shall quantize the fields in a curved background. 
The bosonic fields in any potential U{(l)) obey the Klein-Gordon equation which reads 



1 



(13) 



Setting X = o,4> and V{x) = U{x/a), one gets the equation of motion written in a spacially fiat FLRW metric 

X" - Ax = -a^-^ + -Ra\ 
dx 6 



(14) 



where R is the scalar curvature, a is the scale factor and a prime stands for a derivative with respect to the conformal 
time. The friction term in (|13l) coming from the kinetic term has been absorbed in an effective mass term in (|14|l . 
Since fiat spatial sections are considered, the free real scalar fields may be expanded on a plane wave basis 



ry) 



d^fc 



a{7]) J (27r)3/2 



Uk{ri)e'''-^af: + ul{ri)e-'^-^ai 



a{r]) J (277)3/2'^'= 



(15) 
(16) 



where the operators and at are respectively the annihilation and creation operators and satisfy the equal time 
commutations 



S{k - k'), [a^jflfc,] = 0, 



at, at 
fc' k' 



0. 



The canonical quantization imposes a condition on the Wronskian normalisation 

In the following, we will define the Wightman function to be 

Gk{v,v') = Uk{v)Kiv'), 

its expression in a given background depending on the choice of the vacuum. 
In a FLRW background with fiat spatial sections, the Dirac equation reads 



trd^ (a^/V) - ma (a^/^) = 



(17) 

(18) 
(19) 

(20) 



In the massless case, a free spinor in a spacially flat FLRW spacetime is just conformal to a Minkowski spinor. In 
the following, for the sake of simplicity we shall only consider massless fermions. Moreover, we will consider Majorana 
spinors satisfying the conjugation relation 



4, = a^p) 



(21) 



where C is the charge conjugation matrix that satisfies the property C^^ — — C — . Consequently the decompo- 
sition of a non-massive Majorana spinor on a plane wave basis is trivialy related to the Minkowskian decomposition 
and reads 



1 



,3/2 J (2^)3/2^ 2^ ^"l'^ ^"l'' ^k 

j = l 



(22) 



where k stands for ||fc|| and u^^^{k) et v^^^k) are four-components spinors related by u^^^ = C(v^^^) . The fermionic 
operators commute with the bosonic ones and satisfy the anticommutation relations 



{h\,v%} = 5,y6{k-k') , [hiy-,] - , {v\yl] 



0. 



(23) 
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The free vacuum is defined by the conditions 

aj:\0) = bi\0) = V j,k. (24) 
As in Minkowski, the mode functions satisfy the completness relation 

2 2 

Using the relations H21|) and H25|l . one also finds 

4 4 

^ (0|V^J:.')?M^")|0)(0|Va(x')V'/3(x")|0) = - 5] (0|V^J:z:')^M^")|0)(0|V.(:r:')?/3(^")|0) 

a,/3=0 a,/3=0 

1 /d3pid|2^^^_^^(^,^,^„)^.(,-,+,~,).(,'_,") 



a3(r/')a^(?7") 7 (27r) 



with 



Ap, p, V") = !^^^^^^l^e-^(f i+P^^f"'-"") (27) 

4pip2 

where we set pi = \\pi\\. 

We are now ready to calculate radiative mass corrections. In the next section, we set up the formalism we shall use 
to study one loop mass correction in Yukawa potentials and compare our results to the effect of a small mass. 

III. MASS CORRECTION 

A. Quantum level formalism 

In order to study radiative mass corrections in a generic spacially flat FLRW background, we are interested in the 
calculation of the two-point function 

,,M<t'k, imWim , ?72)0fe, (?72) |0),, (28) 

where is a scalar field, |0),,. is the vacuum state at the time rji and ^7(771,772) is an evolution operator that makes a 
state evolve from the conformal time 772 to the conformal time rji . 

In a free field theory, this expression may be easily computed. For example in an Heisenberg picture, the vacuum 
state does not evolve and the form of the field depends only on the background and on the choice of the vacuum state. 
However this expression must be modified in an interacting field theory and the corrections can only be evaluated 
order by order in the coupling constants in a perturbation theory approach. 

In an interaction picture, the states evolve under the action of the interactions while the fields evolve as if they 
were free. In this representation, the vacuum state evolves as 

|0), = (7(77, 770) |0) (29) 

where |0) is a free vacuum at some time 770, assumed far in the past. The two-point function becomes 

,,(O|0£^(77i)t/(77i, 772)0£^ (772)10),,, = (0|C/(77i,77o)-Vfc,(^i)t/('7i,^2)^fc,(7?2)C/(772, 770)10) (30) 

where the operator [/(771, t^q)"^ is just the inverse of the evolution operator C/(r/i,77o), that is to say (7(770, 771). 
One can find an explicit representation the operator U{rji, 772) which reads in the interaction picture 

-i I dri' Hi{ri') 

<J Tin 



U{i]i,r]2) = ^ 

where Hi{rj) is the interaction hamiltonian and 771 is a conformal time greater than 7/2- 
One has to be careful with time ordering since, for 771 < 772, 



C/(?72, 77i) =A 



+t dri'H,(ri') 

g J ?72 



(31) 



(32) 



where A stands for the anti-time ordering operator. 
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B. Effect of a small mass 

Let us study the correction to the non-massive scalar mode function from a potential of the form 



(33) 



where Sm? is assumed to be a small parameter. 

As a first approximation, we can decompose the scalar field on a plane wave basis, taking into account the loop 
correction. The field then reads 



(27r)3/2 



where we defined 



(34) 



(35) 



The mode function uf^ is the free mode function whereas is the first order correction in Sm^. 

Expanding the mode function in power of the small parameter 6m^ in the equation of motion, one gets the system 



,"(0) 



fc^ ]u 



,(0) 



that could be solved using the retarded Green function 



(36) 
(37) 

(38) 



where Q is the Heaviside function. In the last expression, we used the normalisation p8f) . To first order in dm'^, the 
correction to the mode function reads 



Expanding the Wightman function as 



and inserting H39|l into H4U|) . one then obtains 



,(0) 



,(1 



(0)' 



,(1)' 



(39) 



(40) 
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+ 2 / di^'5m'a'{T^')Gf\m,^')l^ Gl'\m,rj') 

JV2 



(41) 



This correction to the free propagator should be compared to the divergences which may be induced by other 
potentials such as Yukawa couplings of the Wess-Zumino model. 



C. Radiative mass correction from a bosonic loop 

Let us now compute the one loop scalar self-mass induced by a quartic potential 

A 



with the related interaction Hamiltonian 



H,{r^) = j d'x^g^ct>\x) = j d'x-^ 



(42) 



(43) 



Expanding the evolution operator to first order in A, we get from eq. (|30|) 

nMxuM)XkM^\0)v. ^ - » / d77'(0| (r/2), |0) 

•'rjo 

- I rdv'{0\ [x*:, ('71), i?/(^')lx£, ('72)10). 

•'V2 

Using the Wick theorem, this expression transforms into 

= 6A / drj' / d^f'Im (0|Xfc^(r/i)x(a^')|0)(0|Xfc,('72)x(a;')|0)(0|x(a;')x(a:')|0) 

J 1)0 J 

+ 6A rdv' I d3f'Im[(0|Xfe^(r/i)x(a:')|0>l (0|x(x')Xfc,(r/2)|0)(0|x(x')x(a^')|0). 
Finally, using the definition p9|) one finds the contribution from the bosonic loop 

m{0\XkSm)Xk^mmm = 6A<5(fci + fc2)^''dry'Im[Gfc(r;i,77')Gfe (7,2,77')] J -^Gp^^v') 

Comparing eas. H41() and H46|l . one finds a mass contribution from the quartic interaction at one loop order 



(44) 



(45) 



3A r d^p 

This result holds in any spacially fiat FLRW spacetime. 



(46) 



(47) 



D. Scalar self- mass from the fermionic loops 



We now pay attention to the contribution from fermionic loops through a Yukawa coupling of the form 

where ip is a four component Majorana spinor. The interaction Hamiltonian is given by 

A 



Hiiv) = / d'xV^^q^{x)i>ix)^ix). 



(48) 



(49) 



The two-point function (|30|) involves three different evolution operators namely U{'r]2, 770), U^^ivi^ Vo) ^^'^ U (ryi, 772). 
The first non vanishing contribution - after the tree level - comes from the expansion of those operators to second 
order in the coupling constant. Expanding the first two as 



t/(?72,77o) = l-« / dri'Hiii)- drj' dr]" Hiir]')Hiirj") + 0{X^) 

U-Hrii.m) = ^ + ^ d77'iJi(77') - / d77' / dr7"i/i(r7")fl"i(r?') + 0(A3) 
leads to three different terms in the two-point function (|30|1 

V2 rv' 

dr;' / dr^''{0\cj,j:^im)(bj:^im)Hi{v')HiW'm , 

10 J "no 

"'dr;' r dr,'\mi{i')HM)hS'^i)hM)\^) ' 

VI fV2 

d^' / d,7"(O|i?z(77')0£,(m)0,-,(^2)i?/(r7")|O) 

vo •'■no 



(50) 
(51) 

(52) 

(53) 
(54) 
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whereas including the expansion of U{rji, 7)2) exhibits three other terms 

dyy' / d7;"(O|0,,^(77i)i?/(r7')i?/(^")'^fc,(^2)|O), (55) 

ri2 Jri2 

m rm 

dry'/ d,7"(O|i//(r7')0£^(??i)i//(?7")0fc,(^/2)|O) , (56) 
dry' / d,7"(O|0,,^(ryi)i?/(ry')0£^(ry2)i//(r/')|O). (57) 

')2 »?0 

All those terms describe fermionic loops and involve contractions over spinorial operators as well as integrations 
over running momenta. Such expressions - introduced in section II (see eqs.(j2nj and (j2Z|l) - only depend on the 
following quantity 

/(r7',ry"; A,fc) = j ^ j dV2Ap„p. (ry', r?")<5(pi + - k), (58) 

where Apj.p2 (ry', ry") is given by ea. H27() and A is a ultraviolet three-dimensional cut-off. The integration in the 
infrared sector does not lead to any divergence as it could be expected from the spinors being conformal to Minkowski. 
Summing all the terms (|52I57(I and noting that /*(?]', ry"; A, k) = f{ri", rj'; A, k), the two-point function may be written 
in a simpler form after cumbersome calculations 

m(0|Xfe,(^i)Xfc,(^2)|0),, = iX^S^h+k2)l 



dry' / dry"{lm[Gfc(ryi,ry')]Im[Gfe(ry2,ry")/(ry',ry";A,fci)] 
+ Im [Gfc(772, ry')] Im [Gk{vi,v")f{v\ v"; A, ^i)] } 
+ /"dry' /' dry"G^(r72,7?")Ini[Gfc(r7i,7y')]Im[/(ry',ry";A,A;))] 



rm rn^ \ 

+ dry'/ dry"Im[Gfc(ryi,ry')]Im[Gfe(ry2,ry")]/(ry',ry";A,fc) . (59) 

Jn2 Jrio ) 

As low running momenta do not contribute to /(?]', ry"; A, fc), there is no infrared contribution to the bosonic masses. 
Thus the radiative corrections come from the ultraviolet regime, i.e. from large momenta. The expression (EH being 
a product of a divergent function by an oscillating one, the main contribution occurs when the conformal times are 
close to each other. We can then expand the expression (|59|l with respect to the infinitely small parameter e = 77" — 77' 
and the expansion of the Wightman function reads 

Gfc(ry, ry") = Gfc(ry, ry') + (r^' - ry')GU^, ry') + ^^^^^^Gl'(ry, ry') + O [(ry" - ry')^] (60) 
where ' stands for a derivation with respect to the second variable. We also define the real quantities indexed by j 

I,iA,k) = i-^y-'l' dry"i^^^-y^ J§^J dV2Ap„,,(r7',ry")5(pl+p2-fci) 



dV2 



(2^)3 y "^(pi+p2)^+i 



^ _ Pi ■P2 
P1P2 



SiPl +P2-k) (61) 



to get the expansion 



Vo 



dry"Gfe(ry,ry")/(ry',ry";A,A:) =5]z^"-iGi^')(ry,ry')/,(A,fc). (62) 

j 

The computations of Ij{A, k) for j = 0, 1, 2 are sufficient since Ij{A, k) converges for j > 2 and give 

A2 fc2 /A\ ^^/k 



^o(A,fc) = ^-^^^(t)+o(^) , (63) 
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1 



A 



Using the equation of motion for the Wightman function in the background metric 

a 

and inserting ea. H62|l into ea. l|59() . the loop correction reads 

A2 



(64) 
(65) 

(66) 



m(0|Xfc,(^i)Xfe,(^2)|0>^, = -^<5(fci+fc2){ 



/ dr,' 
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drj' 



_ — In 



2a 



G 



— Re[Gfc(77i, 772)]ln - 



47r2 



AX' 



A 



'?2 



dry' / dv"lm[Gk{vi,v')n^[Gk{v2,v")]f{v\v"-:^,k). 



(67) 



Let us compute the last term. Its main contribution comes from rj" and rj' being close to 772 . We can then expand the 
Green function with ?]' — 772 ^ and 77" — 772 ^ 



Im [Gfc (771 , 77')] Im [Gfc (772 , 77")] = E [^fc^ ('^i ' ^2 ) 

J,; 

We are left with two terms of the form 



Im 



G'i\m,m) 



iv' - mV iv" - mY 
7! II 



no 



(Pl +P2)J + 1 ' 

(pi + P2y+^ 



-i(pi+P2}iV2-Vl) 



(68) 

(69) 
(70) 



The equation H70() provides us with a term whose phase is proportional to the difference of the conformal times. Since 
this term oscillates rapidly when |?7i — 7^2! > 1/A, we can neglect its ultraviolet contribution unless the times are equal, 
which would lead to a trivial result for H7U|) . 
The last term of eq. (|67f) then becomes 



r,2 



d77' / d77"Im [Gfc(77i, 77')] Im [Gfc(772, 7;")] 7(77', 77"; A, k) 
i?2 J no 

= ^Im [g*^'^ (771, 7,2)] Im [gI'^ (772, 772)] (z)'--'/,+,+i(A,fc). 

The only divergent term is obtained when / = 1 and j — and reads 

- i Im [Gk im , ?72)] Im [G'^ (7/2 , 7^2 )] /2 ( A, fc) = - ^Im [Gfc (771 , 772 )] /2 ( A, /c) 

where we used eq. (fTH|) . 
Finally, we are left with 



(71) 



(72) 



m(0|Xfc,(?]i)Xfc,(%)|0),„ = 2^{ki+k2)< 
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Vo 

+ / drj' 

J 1)2 



A2- 

a2 



a" 


/ A 




















U 











Im[Gfc, (771, 77')Gfc, (772,7,')] 
GI.^(772,77')Ini[Gfc,(77i,77')] 



(73) 



There is no linear divergence left, in agreement with the local Lorentz invariance. The additional term may be 
cancelled by a field strength renormalisation setting 



with 



Vk{ri) ^ (1 + 5Z)vk{ri). 



(74) 
(75) 



where /i is a typical energy scale depending on the renormalisation condition. /^From ea. l41() . one can see that the 
fermionic loop gives a mass to the scalar field 



X2 











2a ^ 





(76) 



One can remark a background dependent logarithmic divergence proportional to the scalar curvature. 
In the case of a Yukawa coupling l|48|) . the propagator of fermions is just modified by a trace over an odd number 
of gamma matrices which vanishes. Fermions do not acquire a mass through a Yukawa coupling. 

IV. INTERPRETATION 

In the context of the Lagrangian the one loop scalar masses read from eas. (|47|) and H76|l 

A2 



27r2a2(?7) 



4 / dpp2G,(7y,7?)-A2 + ^lnf- 
J 2a \^ 



(77) 



In the following, we compute this expression for different FLRW spacetimes and special attention is paid to de Sitter 
spacetime. 

A. Flat spacetime 

In the flat case, the exact computation of the radiative mass corrections in the massive case may be easily dealt 
with and one expects to recover the Minkowski results by setting a — 1 and 



1 



1 



2uJk 2VWT' 



(78) 



Considering a non- vanishing mass entails an additional correction from the trilinear coupling in the bosonic sector. 
Finally, one finds a zero mass correction for both the scalar fields and a field strength renormalisation as it should be 
expected from a supersymmetric model in Minkowski spacetime [isl Il7j| . 



B. Generic flat FLRW spacetime 

To compute the expression (I77|) in a generic flat FLRW spacetime, we need to specify a vacuum since the choice of a 
vacuum state is not imposed ^] . The simplest choice consists in a Bunch-Davies vacuum whose limit is Minkowskian 
far in the past. With this convention, for any spatially flat FLRW spacetime with constant acceleration, the mode 
function is given by (isj 



Uk{ri) 



irrx' 
"IT 



(79) 
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where 



1 dH 



1 dH 



' aH^ drj dt ' 1-y 



a ■ 



k 



(80) 



and where H^^^ is the Hankel function of the first kind with the index ly. We will only pay attention to period with 
constant acceleration, i.e. y = est. To study the ultraviolet divergences, we shall use the following expansion (see e.g. 
d) 



E 



r(i/ + n+ i) 



Expanding the expression (|79|l for large momenta, one gets 



Uk{ri) 



j^i(ra:— r7r/2 — 7r/2) 



2k 



r+l (r + 2)(r2-l) 



which leads to 



If r + 1 

Gfc(r/, 77) = Uk{7l)ul{ri) = ^ ( 1 + 1^ 



2k 



2a /c2 



(81) 



(82) 



(83) 



Inserting this result into (|77|) , one obtains the one loop mass correction in any constant acceleration FLRW spacetime 
in the ultraviolet regime 



a 



A* 



01 1 



(84) 



where /j, is a typical energy scale. 

There is no quadratic divergence left since we should recover a supersymmetric action on small scales. The local 
Lorentz invariance prevents any linear divergence. The only remaining one is due to the curvature which is in 
agreement with (|11|) and follows directly from a WKB approximation from the equation of motion (|14|l . 

The infrared regime may also lead to divergences. We will only pay attention to the infrared divergence due to the 
bosonic loops since low momenta from the fermionic loops have no contribution. Using the following expansion in the 
case of a non-integer index 



r(i + z.) 



(I)" 



1 



sin(7riy) r(l - I/) V2 



the Wightmann function reads in the small momentum limit 



AaH 



(r/2)^'^ k •,2a+r-l+2v 



+ 



1 



(r/2) fc -.2a+r-l-2iy 



r(l -iy)'aH 



(85) 



(86) 



Inserting this expansion into eq. H47|l gives a vanishing contribution if both the integrals converge, i.e. if 

2a + r + 1 + 2i^ > -1 and 2a + r + 1 - 2i/ > -1. 

Hence, in the massless case, there are infrared divergences for any > 3/2. For a scale factor growing as a{t) ~ t^^, 
we get from ea. l|5n|) 



3/3-1 
2(3-2' 



(87) 



For power law expansion such as /? > 2/3, the scalar mass gets an infrared divergent contribution. In the limit of a 
massless scalar field in a de Sitter space, v exactly equals 3/2 and an infrared logarithmic divergence does appear. 
This infrared behavior has been investigated in 
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C. Consequences on de Sitter inflation 



In this subsection, we apply the previous results to the specific case of de Sitter spacetime which is assumed to 
properly describe the phase of inflation in the early universe. In particular, light scalar fields are bound to generate 
significant fluctuations during inflation if their masses are small compared to the Hubble scale. We want to check 
that their masses are protected against one loop mass corrections. 

In the following, we shall consider a coupling to the curvature which may be added to the Lagrangian (^). The 
curvature in a de Sitter spacetime is constant and given by R = 12H^. Hence, the additional potential looks like a 
mass term for the scalars 

Using the previous results, the effective masses found in (jllll read 



'off 



(89) 



This effective mass is constant and vanishes in the case of conformally coupled scalar fields. To express the quantum 
corrections to the masses, let us choose the Bunch-Davies vacuum for which the mode function reads 



(90) 



where H^p is the Hankel function of the first kind and — 9/4 — 12^. Applying the previous expansions of Hankel 
function, the self mass in the ultraviolet regime reads 



^.^^ = ^(1-60 In 



(91) 



Not surprisingly, there is no additional ultraviolet logarithmic divergence for a conformally coupled scalar field 
with C = 1/6. 

The de Sitter solutions for the massive mode functions are well known and read in the Bunch-Davies vacuum 



Ukiv) 



(92) 



with i''^ = 9/4 — m? /H^ — 12^. From eqs. 147|) and (|85|l . we can see that the contribution to the mass from the 
infrared regime vanishes for any index such as |Re(i^')| < 3/2, i.e. as soon as m? + 12^H^ > 0. 

These results may be applied to the minimally coupled case for which the Wigthman function takes the simpler 
form 



1 



(93) 



The radiative mass corrections read 



— 



Snip = 



27r2c 



27r2a2 



+ 2a^H^ In 



A" - a'H' In ( - 



A 

Air 
A 



Sm — Suit, + 5m-p = 



27r2 



In 



■In 



Air 



(94) 
(95) 
(96) 



For the light scalar fields to develop a significant fluctuations one should get rid of quadratic divergences - as it is 
the case in supersymmetric models - and impose a condition to the coupling constant such as the mass is protected 
against radiative correction 



A^lnl - 1 «1 



(97) 
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which could be easily satisfied for a coupling constant smaller than 10~^. 

The infrared divergence that appears in the mass corrections is related to the l arg e scale structure of spacetime 
and have been extensively discussed in the literature H, lH, 123, l2l, 12^ l2y, 123 ■ Such a divergence naturally 

appears in the real space propagator, the cut-off being of order of the initial comoving horizon, i.e. : 

Air ~ aoH, (98) 

where ao is the scale factor at some initial time [s^- In the context of de Sitter spacetime, the transition between a 
quantum and a classical behavior occurs at the horizon crossing : perturbations modes leave the horizon and behave 
classically outside the horizon 0,11^ 113 ■ 1^ then natural to pay attention to a classical description of the divergences 
which is the aim of the following subsection. 

D. Connection to a classical viewpoint 

During inflation, the universe exponentially grows. The equation of evolution for a stochastic classical field which 
describes the large scale perturbations is given by ea. (|13|l in the low momentum limit 

6s" + 2H5s = -^a^S{5s) (99) 

where S stands for the source term. In the case of a parabolic potential (|33|l with a quartic self interaction H42() . it 
takes the form 

S{5s) = 5m^5s + \5s^. (100) 
In a pertubative approach, we can expand the field into powers of the coupling constant 

5s = 5s(")+,5s(i) + ... (101) 
and define the Fourier transform of the field for each order of perturbation 

As a stochastic gaussian field, has the following properties 



(0)* (0) 
k -k 

{ufuf) = 6^{k + k')P{k) (104) 



u'y> = u'"L (103) 

k —k 



where the brackets stand for an ensemble average and the power spectrum, P{k), entirely defines the field distribution. 
To zeroth order, the modes are frozen and we get from the equation of motion 

= est. (105) 

The first order is sourced by the free field 

<5s(i)" + 2H5s'^^y = -a25(5s(0)) (106) 
where all the time dependence in the r.h.s. is expressed in the factor a? . We then get 

5s(^) (.f,77i) = -S'(^s(°))e~^ / d772 a(?72)^e^. (107) 



1)0 



In the case of a parabolic potential ()33|l . the first order mode function becomes in momentum space 

uf{r]) = -5m^I{'q)uf (108) 

where we defined 

/(?/)= / d77ie "("i) / d??' a(77')^e (109) 

J 110 Jrio 
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^ifc('?i)^ifc,(?72)) -<53(fc + fc')P(fc)- / ^P{p). (110) 



Inserting eas. H103|l . (|104|l and H108|l into (gOl, one gets 

{2n) 

In the case of a quartic potential (|42|l . one finds in momentum space 

and the two-point function becomes 

K^(,7i)u,,,(7y2))=<53(fc + fc')P(fc)-3A<53(fc + fc')m[/(^i)+/(^2)] / 7??T^b)- (112) 



(27r) 



The correction to the mass from the quartic interaction reads 

d^fc 



Sm' = 5Xj ^P{k) (113) 
which is the same result as in the quantum level calculation since 

P{k) ^ ^G,irj,v) (114) 

with k going to zero. In the case of a Harrison-Zeldovitch spectrum, i.e. P(k) ^ k~^, one recovers the same infrared 
logarithmic divergence as in the previous subsection. 

V. CONCLUSION 

We calculated the radiative corrections to the masses of test scalar and fermionic fields in a Wess-Zumino toy model 
embedded in a spatially flat FLRW spacetime. Our results emphasized that the curvature acts as an effective time 
dependent mass m^g = —R/6 for bosons whereas fermions stay massless. In a de Sitter spacetime, this effective mass 
is constant and given by the Hubble scale. 

The final results show that radiative corrections exhibit divergent contributions to the scalar masses both in the 
ultraviolet and infrared regimes. The leading ultraviolet divergences, quadratic and linear in the energy cutoff, 
cancelled both as expected since one should recover the Minkowskian structure in any locally inertial frame. However 
a subleading ultraviolet divergence, logarithmic in the energy cutoff and proportional to the spacetime curvature, 
survives. If it gives masses to scalar fields, it docs not however render them too heavy compared to the Hubble scale. 
Hence, no significant fine tuning will be required for those fields to develop superhorizon fluctuations. Moreover, the 
infrared divergences that we observe arise as soon as the scale factor grows faster than t^/'' . They can be interpreted 
as those due to a self-interacting stochastic classical field. 

These results call for some comments. 

First, one should carefully deal with supersymmetric models in an expanding universe. In supergravity theories, 
the behavior of the fields - including the graviton and the gravitino - has to be deduced from an action which is 
locally supersymmetric. This task is often difficult to handle with. Imposing an expansion of the background, i.e. a 
dynamics for the graviton, and setting the vacuum expectation values of the fermions to zero would generally lead to 
a supersymmetry breaking as confirmed by our result. 

In this context, the use of the Coleman- Weinberg formula in case of an expanding universe can be questioned. 
For instance, in inflationary models - such as F or D-term inflation models - the one loop effective terms due to the 
curvature should be compared to the usual one loop terms of the flat spacetime limit jSOj. Although one does not 
expect this contribution to deeply affect the dynamics at one loop order, consequences of such terms should probably 
be studied in more details (3lll3^ . 

Finally what we have shown here is that light test scalar fields do not acquire a contribution larger than the 
Hubble constant from one loop corrections when they are embedded in chiral multiplet. This result extents what 
was known for the inflaton field. As a consequence it is reasonable to assume that light scalar fields more or less 
generically develop superhorizon isocurvature fluctuations. And then, as was shown in 0, nothing prevents those 
from containing significant non-Gaussianities although whether it is generically the case is still an open question. It 
anyway puts non-Gaussian inflationary models on a much firmer ground. 
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